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We discuss a new mode of ionization front passage in semiconductor structures. The front of 
avalanche ionization propagates into an intrinsic semiconductor with a constant electric field 
in presence of a small concentration of free nonequilibrium carriers - so called preionization. We 
show that if the profile of these initial carriers decays in the direction of the front propagation with 
a characteristic exponent A, the front velocity is determined by ii/ ~ 2/3m/A, where /3m = /3(i?m) 
is the corresponding ionization frequency. By a proper choice of the preionization profile one can 
achieve front velocities Vf that exceed the saturated drift velocity Vs by several orders of magnitude 
even in moderate electric fields. Our propagation mechanism differs from the one for well-known 
TRAPATT fronts. Finally, we discuss physical reasons for the appearance of preionization profiles 
with slow spatial decay. 

PACS numbers: 85.30.-z,72.20.Ht,71.55.-i 



Propagation of impact ionization fronts in semicon- 
ductor structures represents a spectacular nonlinear 
effecljii^i^ which has important applications in pulse 
power electronicsi^ In reverse-biased p'^-n-n'^ diode 
structures ionizing fronts propagate faster than the satu- 
rated drift velocity Ugjii^i^i^ Such superfast propagation is 
possible due to the presence of small concentrations no,j3o 
of free electrons and holes in the depleted region. These 
free carriers which initiate an avalanche multiplication 
are often coined as "pre-ionization" of the medium. ^'^ Ac- 
cording to the conventional concept of ionization fronts in 
TRAPATT (TRAped Plasma Avalanche Triggered Tran- 
sit) diodes^i^ the avalanche multiplication occurs within 
the ionization zone of length If — eeo(£'m — E\j)/qNd 
where electric field exceeds the effective threshold of im- 
pact ionization Eb (Fig. 1, curve 1). This length is fi- 
nite due to the slope of the electric field in the n base 
dE/dx = qNd/eeo which depends on the doping level Nd 
(note that no,po ^ N^)- The finiteness of the ionization 
zone if prevents a uniform avalanche multiplication in 
the whole n base and thus ensures the existence of the 
traveling front mode of avalanche breakdown. However, 
this concept is not applicable to p-i-n structures with in- 
trinsic {Nd = 0) base (Fig. 1, curve 2) as well as to short 
overvoltaged structures because in both cases E > Eh 
in the whole n base (Fig. 1, curve 3). On the basis of 
TRAPATT-like front concept one would expect that in 
these two cases pre-ionization of the high-field region trig- 
gers quasiuniform breakdown ruining the traveling front 
mode. 

In this paper we argue that superfast impact ionization 
fronts are nevertheless possible in p-i-n structures where 
E > Eh everywhere in the high-field region providing the 
concentration profile of initial carriers no (x) , po (x) decays 
in the direction of front propagation. The propagation 
mechanism of such front is completely different from the 
conventional TRAPATT-like front. We find the front ve- 
locity analytically and show that it is controlled by the 
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FIG. 1: Electric field profiles E{x) in the traveling ionization 
front. Profile 1 corresponds to the conventional TRAPATT- 
like front with a finite size of impact ionization zone £ / . Profile 
2 corresponds to p-i-n structure (Nd — 0). Profile 3 corre- 
sponds to p'^-n-n'^ structure with low n base doping. 

slope of pre-ionization profile, and that it can exceed Vg 
by several orders of magnitude. 

We consider a planar impact ionization front and de- 
scribe it by the standard drift-diffusion mode%S. which 
consists of continuity equations for electron and holes 
concentrations n and p and the Poisson equation for the 
electric field E. For a self-similar front motion with con- 
stant velocity Vf these equations can be simplified by 
introducing new variables a = n-\-p, p = p — n.^ For 
intrinsic semiconductor {Nd — 0) the equations for cr, p 
and E in the co- moving frame z — x + Vf t become 

^[vfC7 + viE)p]-D^^2aiE)v{E)a, (1) 

v{E)c7 + Vf p- D^^O, (2) 

dE q 
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FIG. 2: Concentration of electron-hole plasma api, normal- 
ized by (TO = 2eeoaoEo/q, behind the front as a function of 
electric field E^n for Townsend's approximation of the impact 
ionization coefficient a{E) = uq exp{~ Eo / E) . 



where v{E) is the drift velocity, D is the diffusion coeffi- 
cient and a{E) is the impact ionization coefficient. Here 
we neglect recombination and assume that electrons and 
holes are identical in a sense that a{E) = an{E) = ap{E) 
and v{E) — Vn{E) — Vp{E). For an infinite domain the 
boundary conditions are E — > E^, tr, p — > for z —^ ~oo 
and E, p 0, a CTm for z +oo. 

In the simplified case of D = we use Eq. ^ to ex- 
clude p from Eqs.([T|) and Q. This yields 
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= 2 a{E)v{E)c7, 
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Then by dividing equations ([4]) and ([5]) and integrating 
over E we immediately find the dependence <t{E) in the 
moving front: 
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t{E) dE. (7) 



q v)-v^{E)jE 
The plasma concentration far behind the ionization zone 
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depends only on the electric field £'m (Fig. 2). The de- 
pendences p{E) and n{E) in the traveling front 
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q Vf ± v{E) 



a[E)dE. (9) 



follow directly from Eqs. (|6l7p and are shown in Fig. 3. 
Remarkably, a traveling front solution exists for any Vf> 
Vg . Within the approximation D = the slowest solution 
corresponds to the shock front (discontinious at E ^ Ey^) 
that travels with a saturated drift velocity Vf = v^. 




FIG. 3; Concentrations of electrons n and holes p, normal- 
ized by (JO — 2esoaoEo/q, in the traveling front as func- 
tions of electric field E for different front velocites Vf/va = 
1,1.1,1.2,5,10 (curves 1,2,3,4, and 5, respectively) and ap- 
proximations a{E) = ao exp{— E / Eo) , v{E) = Vs/{E + Eg); 
Es/Eo = 0.01. The applied electric field is E^ = Eo/2 = 
50Es. Note that for the chosen v{E) approximation v{Em) = 
0.98 lis, hence vf > Vs for all shown curves. 



The above analysis does not allow to select a physi- 
cally relevant solution and hence to find the actual front 
velocity Vf. The selection problem remains in the case 
oi D ^ 0. This is a general feature of fronts propagat- 
ing into linearly unstable state [see Ref. and references 
therein]. Ionizing fronts belong to this class since the 
state {E = (T = 0) is unstable: due to i^m > Ei, any 
amount of free carriers leads to avalanche multiplication. 
It has also been suggested and confirmed by numerical 
simulations that in gasesiS^iil and semiconductors^^ ion- 
izing fronts are so called pulled fronts. For a pulled front, 
the dynamics in the part of the front where avalanche 
multiplication and screening are essentially nonlinear is 
subordinated to the linear dynamics of the front tip which 
fully determines the propagation velocity.— The dynam- 
ics of the front tip is described by the linearized (near 
the state E = E^n, cr = 0) version of equations I|ll2p with 
constant coefficients v{E) = Vs and a{E) = a{Em) = am 
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FIG. 4: The dispersion relation Vf{X). 



FIG. 5; Linear marginal velocity and the steepness A* 



Here we take into account that in ionizing fields v{E) = 
Vs — const. Solutions of these linear equations are expo- 
nential functions (t{z), p(z) ~ exp(Az), where the disper- 
sion relation Vf{X) remains to be found. 

It is known from the theory of pulled fronts that the 
actual front velocity strongly depends on the type of ini- 
tial conditions)^ All possible initial conditions split in two 
classes that lead to qualitatively different dynamics. Lo- 
calized conditions correspond to initial profiles a{x,t = 0) 
that are steeper than the profile exp(A*a::) with a cer- 
tain characteristic exponent A"^: cr{x) < Cexp(A*a;) for 
X ^ — oo, where C is an arbitrary constant. In this 
case the front profile eventually becomes smoother and 
asymptotically reaches the profile cr(z) ~ exp(A*z) at 
z — > — oo that propagates with linear marginal stability 
velocity v* = Vf{X*)^ Any initial profile that is strictly 
equal to zero for sufficiently small value of x also rep- 
resents a localized initial condition. Nonlocalized initial 
conditions correspond to profiles a{x,t = 0) with slow 
spatial decay that do not meet the above mentioned con- 
dition cr{x) < C exp(A'^a;) for x —^ —oo and hence are 
smoother than exp(A'^a;). In this case the front velocity 
is fully determined by a{x,t — 0): for the initial profile 
exp(Aoa;) with Aq < A* the front velocity is given by the 
dispersion relation vq — VflXo). 

The dispersion relation f/(A) follows from the charac- 
terictic equation of Eqs. and pT|) 
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and is explicitly given by (see Fig. 4) 
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The critical steepness A* and the velocity v* correspond 
to the minimum point^ and are given by 




(14) 



+ iA + a^£)A, (15) 



where A = ^{a^lYj^ + a,„€. 



The right branch of the w/(A) (A > A*) corresponds 
to fronts whose velocity increases with steepness A due 
to diffusion. According to the concept of localized initial 
conditions, % these fronts are unstable: their steep pro- 
files eventually relax to the profile with exponential tip 
exp[A'^z] that propagates with the velocity v* . The left 
branch of the Vj{X) dependence (A < A*) corresponds 
to stable fronts whose velocity decreases with A. These 
fronts correspond to nonlocalized initial conditions which 
are in the focus of our interest. 

Charateristic values of the dimensionless parameter 
a-a^l = a-mDjvs are 0.1 for Si and 1 for GaAs devices. 
As follows from Fig. 5, A"^ > am in the relevant interval. 
Physically it means that the front propagating with linear 
marginal stability velocity v* is so steep that the valid- 
ity of the drift-diffusion approximation is questionable. 
This problem disappears for much smoother fronts that 
correspond to nonlocalized initial conditions (left branch 
in Fig. 4) if A < A*. 

Eq. leads to a simple relation Vf/vg — 2ani/X for 
the ionizing front velocity in case of prcionization with 
decay exponent A < A*. For such fast fronts the effect of 
diffusion is negligible; in particular, Eq. (|4I5I6I7I8I9D are 
fully applicable. Although Vf increase with am and hence 
with electric field i^m, it is the ratio am/ A which actu- 
ally counts. This ratio can be arbitrarily large resulting 
in front velocities that exceed the saturated drift veloc- 
ity by many orders of magnitude. It means that a proper 
choice of slowly decaying preionization profile gives the 
possibility to achieve fast front propagation in even mod- 
erate (with respect to Eb) electric fields. However, the 
concentration of electron-hole plasma generated by the 
front passage increases with (Fig. 2). Generally, the 
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electromagnetic limitation Vf < c, where c is the velocity 
of light, may be important: for Vf comparable to c the 
full set of Maxwell equations shall substitute the Pois- 
son equation in the model because the feedback from a 
nonstationary electromagnetic field created by the front 
passage on the front dynamics becomes essential. 

Pre-ionization profiles with slow spatial decay can 
appear due to photoionization by photons from dense 
electron-hole plasma behind the front. In this case, 
A""'^ can be roughly identified as the light absorption 
length. This mechanism is efficient in direct-band ma- 
terials and can be relevant for planar fronts in GaAs 
diode structures^ as well as finger-like streamers in direct- 
band bulk semiconductors!^ Another mechanism is re- 
lated to field-enhanced ionization of deep centers in Si 
p~^-n-n^ structures used in pulse power applications.'^ 
These high-voltage structures possess "hidden" deep lev- 
els - process-induced defects - with low recombination 
activityji^ Preionization of the high-field space charge 



region can be due to field-enhanced ionization of these 
deep centers embedded in the n base.^** This ionization 
is more efficient near the p^-n junction where the electric 
field is stronger. Hence the profile of initial carriers de- 
creases along the n base. The characteristic decay length 
A^^ is expected to be a fraction of the n base width 
W ~ 100 /im. At the same time for low doped n base the 
electric field can be above the ionization threshold ev- 
erywhere. This may result in front velocities that exceed 
Vg by several orders of magnitude. Numerical simulations 
of such triggering process will be presented elsewhere. 
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